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Abstract 

A unified theoretical description of ballistic and diffusive carrier transport in parallel-plane semi- 
conductor structures is developed within the semiclassical model. The approach is based on the 
introduction of a thermo-ballistic current consisting of carriers which move ballistically in the elec- 
tric field provided by the band edge potential, and are thermalized at certain randomly distributed 
equilibration points by coupling to the background of impurity atoms and carriers in equilibrium. 
The sum of the thermo-ballistic and background currents is conserved, and is identified with the 
physical current. The current- voltage characteristic for nondegenerate systems and the zero-bias 
conductance for degenerate systems are expressed in terms of a reduced resistance. For arbitrary 
mean free path and arbitrary shape of the band edge potential profile, this quantity is determined 
from the solution of an integral equation, which also provides the quasi-Fermi level and the thermo- 
ballistic current. To illustrate the formalism, a number of simple examples are considered explicitly. 
The present work is compared with previous attempts towards a unified description of ballistic and 
diffusive transport. 

PACS numbers: 05.60.Cd, 72.10.-d, 72.20.-i 
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I. INTRODUCTION 



Within the semiclassical model, charge carrier transport in semiconductor structures is 
commonly described in two limiting situations. In the ballistic limit (carrier mean free path 
much larger than the characteristic dimensions of the sample) , the carriers move unscattered 
all through the sample. The conduction is treated as transmission across the band edge 
profile of carriers emitted and absorbed at external contacts which are represented by large 
reservoirs in thermal equilibriu m 1 ! 2 ' 3 ! 4 ' 5 . This is the picture of ballistic motion thermally 
activated at the contacts in terms of which one derives a current-voltage characteristic for 
the current flowing in the sample. In the other, the diffusive, limit (mean free path much 
smaller than the characteristic dimensions of the sample), one has the Drude model of 
carriers moving from one equilibrium state to another, infinitesimally close-lying one, each 
characterized by the local quasi-Fermi level and by the electric field provided by the band 
edge potentia l 6 ! 7 ' 8 ! 9 ' 10 ' 11 . This yields the local drift-diffusion equations for the current density 
in the sample. 

In the present work, we consider the general case where the ratio between mean free 
path and the characteristic dimensions of the sample may assume any given value. For 
definiteness, we treat electron transport; the case of holes is analogous. We develop a uni- 
fied description of ballistic and diffusive electron transport in parallel-plane semiconductor 
structures, in which the ballistic motion is thermally activated neither only at the distant 
external contacts nor only at infinitesimally close-lying points, but at separations anywhere 
between these two limits. No restrictions are imposed on the shape of the conduction band 
edge profile. 

The unified description relies on the introduction of a "thermo-ballistic current", where 
the electrons move ballistically in the band edge potential over the length of certain randomly 
distributed "ballistic intervals" , at the ends of which they are thermalized. The probability 
for ballistic (unscattered) motion between any two end points ( " points of equilibration" ) 
is governed by the mean free path. The thermalization is effected by the coupling of the 
electron motion to the background of impurity atoms and electrons in equilibrium. 

The equilibration at the end points of the ballistic intervals, i.e., at the points separating 
adjacent intervals, is viewed as absorption of the electrons impinging from either side, and 
re-emission with a momentum distribution determined by the local quasi-Fermi level of the 
background. The absorbed and re-emitted currents do not, in general, counter-balance, 
and produce a background current whose local increment is equal to the negative of the 
increment of the thermo-ballistic current. Therefore, the sum of the thermo-ballistic and 
background currents is conserved, and is to be identified with the physical current. From 
the condition that the integral of the background current over the length of the sample 
vanishes, an integral equation is derived for the quasi-Fermi level inside the sample in terms 
of a "resistance function" x, which yields the thermo-ballistic current and whose value 
X at one end of the sample determines the current-voltage characteristic. This quantity 
depends on the mean free path I and on the shape of the conduction band edge profile E c , 
X = xi}i [-£<;])> an d is called the "reduced resistance". It is the central quantity of the unified 
description of ballistic and diffusive electron transport. The derivation is carried out for the 
nondegenerate case; it turns out that essentially the same reduced resistance determines the 
zero-bias conductance in the degenerate case. 

In earlier paper a 12 i 13 i 14 i 15 i 16 , attempts have been made to treat the synthesis of ballistic 
and diffusive transport for particular band edge profiles. A heuristic unified model yielding 
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a closed expression for the reduced resistance in terms of the shape of an arbitrary band edge 
profile has been presented in Ref.— . The present work is directed towards a comprehensive 
description of electron transport, with emphasis on a consistent theoretical formulation. No 
attempt is made here to apply our approach in the analysis of experimental data. 

In the next section, we introduce the two basic elements of our description, the points of 
local equilibration and the ballistic intervals, and derive the probabilities of occurrence of 
the latter. In Sec. Ill, we construct the thermo-ballistic current and establish its connection 
with the physical electron current. In Sec. IV, the integral equation for the function \ is 
derived and the current- voltage characteristic is expressed in terms of the reduced resistance 
X- The ballistic and diffusive limits are obtained as special cases. A discussion of the degen- 
erate case completes the section. Simple examples illustrating our formalism are considered 
in Sec. V. We examine various types of conduction band edge profile: flat profile (constant 
electrostatic potential), profiles rising (or falling) monotonically, and profiles containing a 
single maximum. We present calculations of the reduced resistance x, of the quasi-Fermi 
level, and of the thermo-ballistic current for one-dimensional and three-dimensional trans- 
port. In Sec. VI, we compare our approach with previous attempts to unify ballistic and 
diffusive transport. Finally, Sec. VII contains a summary of the present work and an outlook. 

II. EQUILIBRATION POINTS AND BALLISTIC INTERVALS 

We consider the transport of electrons in a semiconductor sample of the form of a rect- 
angular slab of thickness S, confined between two plane-parallel contacts (cf. Fig. ^). The 
x-axis is taken normal to the planes of the contacts, and the sample is assumed to be uniform 
in all directions perpendicular to the x-axis. We assume that the scattering at the boundary 
planes parallel to the x-axis is specular so that effectively these planes can be regarded as 
infinitely removed. We are, therefore, dealing with a one- dimensional geometry along the 
x-axis but, of course, the electrons move in three-dimensional space. For a two-dimensional 
electron gas, the latter may be reduced to two dimensions, and for illustrative purposes one 
may even consider electrons in one-dimensional motion (in or against the x-direction) . 

The contacts at x = and x = S are viewed^^ as large reservoirs with Fermi ener- 
gies -Ef(O) and E-p(S), respectively, whose difference determines the voltage bias V. As 
a relaxation mechanism inside the sample, we consider only isotropic impurity scattering 
characterized by a universal mean free path I, which depends on the doping concentration, 
but is independent of the electron velocity. In the present geometry, the spatial dependence 
of all quantities of interest reduces to that on the x-coordinate. For notational convenience, 
we introduce the scaled coordinate £ = x/l, that is we measure all lengths in multiples of 
the mean free path I; the scaled coordinate o = S/l denotes the position of the contact at 
x = S. 

We introduce points of equilibration along the £-axis with the understanding that these are 
actually short intervals (with position £ and length d£) representing the three-dimensional 
sheet enclosed between two (infinitely extended) planes perpendicular to the £-axis and sep- 
arated by the distance <i£. The electrons arriving at such points escape into the background 
and are instantaneously replaced with background electrons emitted isotropically into the 
thermo-ballistic current. The density J e (0 of the particle current emitted at a point £ [for 
the charge current density j e (£), one has to multiply by — e] into the right (left) direction is 
given by 

J e (0 = ±v e n(0, (1) 
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where v e = {2nm* f3)~ 1 / 2 is the emission velocity, and n(£) is the electron density at the 
equilibration point £. For a nondegenerate system (we consider the degenerate case later 
on), we have 

n(Z) = N c e-ft E °®- E *W ; (2) 

where N c = 2(2irm* / (3h 2 ) 3 / 2 is the effective density of states at the conduction band edge, 
E c (£) is the profile of the conduction band edge potential, E F (£) is the quasi- Fermi level 
at the position £ in the interval < £ < a, and (3 = Xjk^T. The contacts at £ = and 
£ = a may be regarded as special points of equilibration characterized by true Fermi energies 
-Ef(O) and E F (a). Including these quantities in the definition of Ep(£), we continue to call 
this function the quasi- Fermi level (now defined for < £ < a). 

Between two points £' and £ of local equilibration, the electrons move ballistically in the 
electric field provided by the conduction band edge. These ballistic intervals have an average 
length of the order of the mean free path. In order to determine the probability of occur- 
rence of a ballistic interval of length |£ — £'|, we briefly digress and introduce the following 
probabilities, which depend on the ratio of the distance between the two equilibration points 
and the mean free path, i.e., on the scaled distance |£ — £']: 

1) — the probability pi(|£ — £'|) that an electron emitted at the point of local equilibration 
£' reaches the point £ without scattering [clearly pi(0) = 1]; 

2) - - the probability p 2 (|^ ~ that an electron emitted at a point £' after having 
been equilibrated there within the interval <i£' reaches the point £ without scattering and, 
analogously, the probability p 2 (|£ — £'|)g?£ that an electron emitted at £' reaches the point £ 
without scattering, and is equilibrated there within the interval g?£; 

3) — the probability p 3 (j^ — £'\)d£'d£ that an electron emitted at a point £' after having 
been equilibrated there within the interval d£' reaches the point £ without scattering, and 
is equilibrated there within the interval g?£. 

The probability functions p n (£) depend on the dimensionality of the electron motion via 
the averaged projection of the mean free path I on the rr-axis (cf. below). 

The probability p 2 (]^ — is the product of the probability Pi(|£ — £'|) for motion 

between £' and £ without scattering times the absolute value of the relative change of this 
"no-scattering probability" (which is the negative of the relative change of the scattering 
probability) within the interval [£,£ + d£], 

p 2 (|£ - £1K = pi(|£ - £1) (- ^rff ) = - *'D ' ( 3 ) 

or 

p 2 (£) = -d Pl {i)/di ■ £>0. (4) 

Similarly, we have 

P3(|£-£lK'd£ = -dp 2 (|£-£W, (5) 

or 

p 3 (0 = -drtO/dt ; e>o. (6) 

For later use, we define more generally 

Pn(0 = -rfPn-i(£)M ; £>0, n = 1,2,3, (7) 

or 

/Vpn(O=Pn-i(0)-p„_i(£); n = l,2,3; (8) 
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here, the so far undefined quantity po(0) is found, with the help of Eq. (7), to be given by 

poo 1 poo A T (J) 

Po(0) = I dftPaCO = y I -*P2(x/l) = Y • ( 9 ) 

The quantity (/) is the effective mean free path along the x-axis; it depends on the dimen- 
sionality of the electron motion via the probability function p2(x/l). 

We now return to the probability of occurrence of the ballistic intervals. The latter have 
been defined as intervals of ballistic motion with equilibration (absorption and emission) at 
the two ends. Three different kinds of ballistic interval are to be distinguished: 

1) In the case of the ballistic interval spanned between the contacts at £ = and £ = a, 
the electrons are emitted and absorbed, respectively, at the contacts. The probability of 
occurrence of this interval is given by Pi(o~). 

2) For the ballistic intervals [0, £] and [£,cx], the equilibration occurs in the infinitesimal 
interval [£,£ + d£], and at the contacts at £ = and £ = a, respectively. These ballistic 
intervals occur with probability £>2(£)<^£ an d Vi(p ~ £)^£j respectively. 

3) Finally, for the ballistic interval [£,£'], the probability of occurrence is given by 

Pa(l£-£'IK'rf£- 

Considering the explicit form of the probability functions p n (£) for the different dimen- 
sionalities, we find for one- dimensional motion 

Pn(0 = e- e ; n = 0,l,2,3. (10) 

For n=l, this is the usual "survival probability" for ballistic motio n 10 ! 11 ; the other cases 
follow with the help of Eqs. (JZJ)-©. In particular, po(0) = 1, which confirms the obvious 
result (I) = I. 

For two-dimensional motion, we have to replace the mean free path I with its projection 
I cos 8 on the x-axis (cf. Fig. Q). Then, changing on the right-hand side of Eq. (jlOj) for n = 1 
the variable £ to £/ cos# and taking the average over the angle 9, we obtain 

/"7r/2 poo f>~^t 

Pi(0=/ d9 cos 9 e-^ cose = dt — , (11) 

mW Jo Ji t 2 V^l 1 ; 

where we have introduced the variable t = 1/ cos 9. More generally, we have from Eqs. Q-© 

poo 

Pn(0=j i V-Vt 2 ~ 1 ; n=0 ' 1 ' 2 ' 3 - (12) 

In particular, 

ps(0 = Xo(£) , (13) 
where i^o(£) is a modified Bessel function^. Furthermore, 

po(0) = |, Pi(0) = l, pa(0) = ~, (14) 

and the effective mean free path comes out as (I) = (tt/4) I. 
In the case of three-dimensional motion, we have 



Pl (£) = 2 d9 sin 9 cos 9 e^ /cosd = 2 f d(cos0) cos9e-^ /cosd 
Jo Jo 



rft^ = 2E 3 (£). (15) 
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Here, E n (£) is the exponential integral of order n;— the factor 2 is a normalization factor, 
since 

f d(cos9)cos9 = - . (16) 
Jo 2 

Using Eqs. (0)-© and the properties of the exponential integral, we find 

p n (0 = 2£ 4 _ n (0 ! n = 0,l,2,3, (17) 

and 

Pn(0) = -^; n = 0,1,2. (18) 
The effective mean free path is (I) = (2/3) I. 

III. THERMO-BALLISTIC CURRENT AND PHYSICAL CURRENT 

The concepts established in the preceding section provide the basis for the derivation of 
the net current at an arbitrary point £ (0 < £ < a) inside the sample of length a (cf. Fig. |5J). 
On the one hand, electrons arrive ballistically from the left, emitted at the left contact at 
£ = or at a point of local equilibration £' (0 < £' < £); after having passed the point £, 
they travel ballistically to the right contact of the sample at £ = a or to points of local 
equilibration at £" (£ < £" < a). The corresponding currents are illustrated by the four 
arrows with heads pointing to the right. On the other hand, electrons arrive ballistically from 
the right, emitted at the contact at £ = a or at a point of local equilibration £" (£ < £" < a), 
after which they move on to the left contact at £ = 0, or to points of local equilibration at 
£' (0 <£'<£), as illustrated by the arrows with heads pointing to the left. 

Inside some ballistic interval [£', £"], the current emitted at £' is determined by the electron 
density n(£'), and the current emitted at £", by n(£ // ) [cf. Eq. (JTJ)]. However, owing to the 
presence of the band edge potential -E c (£), part of the current is reflected, so that we must 
introduce transmission probabilities 

T 5 ,(£', £") = e -/3[W,0-£c(f)] ) Tf//( ^ = e H9[W,{")-fictt")] (19) 

for emission at the end points £' and £", respectively, of the ballistic interval; here, £™(£', £") 
is the maximum of the potential profile within [£', £"] (note that _E™(£', £") is symmetric with 
respect to the interchange of £' and £"). The transmission probabilities (|19|) are the classical 
ones; if tunneling effects are important, the former should be replaced with the corresponding 
WKB expressions (cf. Refs.^2^). In view of Eqs. (JTJ) and (J2J), we now have (0 < £' < £" < a) 

Mf, £") = vMOM?, £") = v e N c e-^'^- E ^ (20) 
for that part of the current across [£', £"] which is emitted on the left at £', and 

Me, n = -vMewe, n = -^-^w^-^o (2 i) 

for the part emitted on the right at £". The net current in the ballistic interval [£', £"] is 
therefore 

J(£',£ / = Me,n + ^(£',£") = v e N c e- 0E ^'*'V E ^ - e^ F «")] , (22) 
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which is antisymmetric under the interchange of £' and £". 

We now find for the total (net) electron current passing the point £ (0 < £ < a) 

J(0 = Pi(a)J(0,a)+ f^'P2^-eV(C» 

Jo 

+ f ^'W£V (o, £") + f Q d? f rf£> 3 (£" - n . (23) 

Here, the first term on the right-hand side represents the net current of electrons passing £ 
while moving ballistically between and a [which occurs with probability pi(a)}, the second 
term refers to the electrons in ballistic motion between any point £' (0 < £' < £) and a 
[probability p^io — £')^£'], etc. In Fig. the four terms in expression are represented 
by the double-arrows labelled I to IV. 

We call the current J(£) the thermo-ballistic current, since it is the sum of weighted 
contributions of ballistic currents emitted and absorbed at the contacts at £ = and £ = a, 
or at the local equilibration ( "thermalization" ) points £' and £" inside the sample. For 
analyzing the properties of this current, we consider its (positive or negative) increment 
over the infinitesimal interval [£,£ + d£], 



dJ(0 = J(£ + rf£) - J(£) = mo/dm 



P 2 (£)J(0,£)+ / V? 3 (£ - O^', 
Jo 



d£ 



+ 



P2(a - £) J(£, a) + f rf£> 3 (£" - £)</(£, £") 



rf£ . (24) 



The increment of the thermo-ballistic current is caused by the net influx supplied to it via 
the emission and absorption (equilibration) processes occurring in the interval [£,£ + c?£], 
as depicted in Fig. El The net contribution of emission and absorption in that interval 
to the thermo-ballistic current in the interval [0,£] (indicated by the double-arrow labelled 
I)is represented by the first term in the first pair of brackets in expression (|24p. the net 
contribution to the thermo-ballistic current in the interval [£', £] (double-arrow II) by the 
second term, and similarly for the two terms in the second pair of brackets (double-arrows 
III and IV). 

In order to compensate the local imbalance between emitted and absorbed currents at 
the equilibration points, we introduce a "background current" J back (£), in such a way that 
the total current J tot (£), defined as the sum of the thermo-ballistic current (|23|) and the 
background current, is conserved, 



J tot (£) = J(£) + J back (£) = J = const. 



(25) 



This current is to be identified with the physical current. Since, by its nature, the background 
current is confined within the sample, it must vanish when averaged over the latter, 



S. J*°*(£) = . 



o a 



Now, averaging Eq. ()25|) over the sample and using Eq. 
"global conservation" for the thermo-ballistic current, 

di 



Jo a 



(26) 

we derive the condition of 

(27) 
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This relation expresses the physical current J in terms of the thermo-ballistic current J(£); 
it will be the basis for setting up a formalism which allows us to obtain the current-voltage 
characteristic, i.e., a relation between the current J and the voltage bias V. 



IV. CURRENT- VOLTAGE CHARACTERISTIC 



A. General form of the current-voltage characteristic 



The thermo-ballistic current «/(£) is expressed by relation ()23|) in terms of the two- variable 
function J(£', £"). In order to obtain an expression containing a function that depends on the 
single variable £ only, we use Eq. (}2"2"|) to write the current J(£', £") in the ballistic interval 
[£',£"] in terms of the ballistic current J(0,£) in the interval [0, £], 



./(£',£") = 7(£',£'V(o,£") -7(£",£'V(o, 



with 



J(0,O = v e N c e-P E ^Y EF(0) ~ e^ F(0 ] 
and the "profile function" 7(£',£") defined as 



7(e',o 



-/3[B-(r,n-^(o,r')] 



(28) 
(29) 

(30) 



The latter function is determined by the band edge profile E c , which actually is a universal 
function of the coordinate x, independent of the mean free path I, i.e., E c = E c (x). Therefore, 
as a function of x' and x", 7 = j(x',x"; [E c ]) is also independent of I. For notational 
convenience, we continue to use the shorthand notation E c (£) and 7(£', £") to represent the 
functions E c (x) = E c (lg) and j(x',x";[E c ]) = 7(l£',l£"',[E c ]). We note that the profile 
E c (x) is indirectly connected with I via the dependence of both E c and / on the doping 
concentration in the sample^. 

With the use of Eq. (|28|). expression (|23j) for the thermo-ballistic current now becomes 



m 
j 



Jo 

f Wo- - en*, o + r de^e - ene, o 

Jo L J? 

r^'L(n + /Vp3(£ , -£")7(r,£ , ) 

7f JO 



x(0 



where we have introduced the (dimensionless) function 

x (0 = jJ(o,o, 



x(0 



(31) 



(32) 



which is defined for < £ < a. If we assume, without loss of generality, that E F (0) > Ep(a), 
the physical current J flows from left to right, J > 0, and the same holds for the ballistic 
current between the contacts, J(0, a) > [cf. Eq. (J2HJ)]. Owing to the effect of scattering, 
the physical current must be smaller than the ballistic current, so that we have 



(33) 
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The quantity x denned as 

X = (34) 

is the central quantity in the current-voltage characteristic following from Eqs. (J29*|) [for 
£ = a) and 

j = -ev c N c e-^U^-e- peV ), (35) 

where j = —eJ is the charge current, E p = E™(0,a) — Ep(0) is the overall barrier height 
along the whole sample, and V = [Ep(0) — Ep(cr)]/e is the voltage bias between the contacts 
at the ends of the sample. For j3eV <C 1, Eq. (J33j) reduces to j = —V/R, where the quantity 
R (resistance times area of the sample) is seen to be proportional to Xi 

R=(e 2 v e N c (3)- 1 e^x- (36) 

We therefore call x the reduced resistance, and the function x{0 the resistance function. We 
will see below that, in general, x is n °t simply given by Eq. (|34j) but must be symmetrized 
with respect to the profile E c . 

Performing the integration on Eq. (j3*Tj) according to Eq. (|27jl. we find 

<r-f(*)x(<r)+ f d£K(<r, x(0 = , (37) 
J o 

with the coefficient function 

/(0 = fri(0 + f d£\£ - £')V2{i - 07(£', (38) 

JO 

and the kernel 

m a = -&*(o + « - rwe - 07«, o + f w - OvM' - ane, a m 

Jo 

(£' < £)■ Equation (|37jl is a condition on the function required by Eq. (|77|. In the 

following, it will be reinterpreted as an integral equation for x(0> allowing us to determine 
the reduced resistance x = xO 7 )- 



B. Determination of the resistance function 

The coefficient function /(£) and the kernel /C(£, £') depend on the mean free path Z and 
the band edge profile E c of the system via the profile function 7(£', £") = 7(££', l£"] [E c ]) [cf. 
Eq. (EOJ)], 

/(£) = /(£;U£cD, (40) 
/C(£,0 = ^(^^V,[^c]). (41) 

Apart from the implicit dependence on the profile contained in / and JC, Eq. f|37jl contains 
the position a = S/l of the right contact explicitly. We now allow this quantity to become 
a variable which may assume any value between o = + and a = S/l. Accordingly, by 
replacing in Eq. (J37|) the quantity a with the variable £, we obtain a Volterra-type integral 
equation for the resistance function x(0> 

t - f(0x(0 + f d?K(t O X (0 = . (42) 
Jo 
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By solving Eq. ()42|) in the range < £ < a, we can determine the function x(0 which 
satisfies condition (p?Tj) . 

As a preliminary to the calculation of the function x(£), we look at the behavior of that 
function for £ — > 0. We use 7(0,0) = 1 in Eqs. (|HH|) and ()39j). and expand the right-hand 
sides of the resulting expressions for /(£) and /C(£, £') [cf. Eqs. (|57jl and (jUS|) below] to first 
order in £ and With the help of Eq. 0, we find /(£) « £ and « p 2 (0)(£ - 2f) 

for f , f ' -> 0, so that Eq. (g2D yields 



X(0 + ) = 1 • 

Therefore, the function x(0 is discontinuous at £ = 0, since by Eqs. and 
Equation (J43|) defines the initial value for solving Eq. (J42j) . It implies 



J(0,0 H 



J 



(43) 
X(0) = 0. 

(44) 



since J ^ 0, we have .Ef(0 + ) 7^ -Ep(O) [cf. Eq. (|29|)]. i.e., the quasi-Fermi level for £ = + 
differs from the Fermi level in the contact at £ = by an amount determined by Eq. (j4*3"|) . 

Equation (j4*2"|) is the key result of the formal development of this paper. It is a lin- 
ear equation as long as the band edge profile E c (£), and consequently the profile function 
7(£', £"), can be considered given, as in the case of zero bias. Otherwise, the profile function 
will depend on the solution x(0 which, via the quasi-Fermi level Ep(£), enters the Poisson 
equation for E c (£). Equation (J42j) then becomes a nonlinear equation which must be solved 
self-consistently together with the Poisson equation. 

Once the function x(0 nas been determined for < £ < er, we can obtain the quasi-Fermi 
level Ep(£) as a function of the parameter J by using Eqs. ()29|) and 



£ F (O = £ F (0) + -ln 



J 



v P N, 



x(£) e 



-f3[E F (0)-E?(0,O] 



(45) 



Eliminating in this expression the current J = —j/e with the help of the current-voltage 
characteristic (J3SJ), we can also write Ep(£) as a function of the parameter V: 



£ F (O = £ F (0) + ^ln 



x(0 
x(*) 



(46) 



The thermo-ballistic current J(£) inside the sample is found from Eq. (|3ip. with the physical 
current J entering as an overall factor only. 



C. Interpretation 

Since the resistance function approaches unity as £ — > + , it is appropriate to write it 
(for £ > 0) in the form 

X (0 = 1 + k(0£ • (47) 

From the examples presented in Sec. V, it emerges that > 0, that k(0) is finite, and 
that in the limit £ — > 00, the function tends toward a finite value. 

The first term on the right-hand side of Eq. 1)470. i.e., the limiting value of x(0 f° r 
£ = 2// — ► + , represents the purely ballistic part of the resistance if this limit is thought 
to have been obtained for I —>■ 00 at arbitrary fixed x. On the other hand, if we let 
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x — > + at arbitrary fixed I, this term is to be interpreted as a contact resistance: from 
Eqs. (129)1 and (|44jl. we find that a finite (but small compared with k^T /e) voltage difference 
V(0 + ) = [Ep(0) — Ep(0 + )}/ e across the interval [0, + ] is associated with a finite conductance 
per unit area 



J 




J(0,0+) 


V(0+) 


= e 


V(0+) 



(48) 



(recall j = —eJ), although apparently there is no hindrance to the current in that interval. 
This conductance is the inverse of the contact resistance located in an ideal, "reflection- 
less" contacl^^i at £ = 0; it is the nondegenerate semiconductor analogue of the contact 
conductance for a three-dimensional degenerate electron conductor at zero temperature^, 
G c = (2e 2 //i)(fc|/47r), where fcf = 2m*[E F (0) - E c (0)]/h 2 . The contact resistance is located 
inside the left contact because the resistance function x(0 is associated, via Eqs. (|2*§|) and 
(132)1 . with a quasi-Fermi level Ep(£) determining the current injected from the right, which 
enters the left contact without reflection^. 

Turning now to the interpretation of the thermo-ballistic current J(£), we restrict our- 
selves, for the sake of transparency, to the special case of constant band edge profile, 
E c (x) = E c (0). Then the current- voltage characteristic reads 



j = -ev e N B e-ft E °<- '*- E *<M-± T (l 



-!3eV\ 



(49) 



[cf. Eq. (jSHJ)], where x( a ) is a universal function of a. The quasi-Fermi level is given by 



£ F (0 = E F (0) + - In 



1 _ x(0 n 



— e 



(50) 



[cf. Eq. 



and the thermo-ballistic current J(£) is obtained from Eq. ()31j) as 



j 



Pi{° - " f ^'p 2 (|£ - ^1) sgn(£ - O X(0 

JO 



(51) 



where use has been made of Eq. (J7j). 

Equation (|51|) may be rewritten in the form 



m 
j 



Pi(Ox(o + ) + £<%' Pl (\z-e\)-^x(0 



Introducing the "internal voltage difference" across the interval [0, £], 

V{£) = [E P (0)-E P {We, 

and the thermo-ballistic charge current = — eJ(£), we now have, using Eqs. 
and dSI|l. 

d 



(52) 



(53) 



•/^ _ Pl(0 .-ball 

JlU X(a) J 



Here, 



(54) 

(55) 
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is the current in the ballistic limit, and 



C&0 = 2^Pi(l£-ri)c(0, (56) 

where c(£) = e/zn(£) is the (local) conductivity, with the mobility /x = 2ev e j3{l). In expression 
(154)) for the thermo-ballistic current at position £, the first term represents the ballistic 
component which enters with weight factor Vi{£) / x{°~) ■ The second term represents the 
diffusive component, driven by the internal voltage gradient —dV(£)/d£, with a nonlocal 
conductivity £(£,£')■ m the ballistic limit, we have = j = j , while in the diffusive 
limit C(£, £') = C(x/l, x' /I) oc 8(x — x'), so that 

j(x) = j = -c(x) 4~V(x) . (57) 

In the Appendix, these limits will be worked out for the case of arbitrary band edge profile. 

The starting point of the development in Sec. IV.A has been to rewrite the current 
</(£',£") in terms of the current J(0, £) containing the quasi- Fermi level E-p(£) at the right 
end of the interval [0, £] extending between the left contact at £ = and the position £. 
However, in the formalism of Sec. Ill no preference is given to either contact. Therefore, we 
could just as well have chosen to write the current «/(£', £") in terms of the current J(£, a) 
containing the quasi- Fermi level at the left end of the interval [£, a] extending between the 
position £ and the right contact at £ = a. The modifications caused by this alternative 
choice are the following. 

The function x(£; [E c ]) is replaced with the function Xi(& [E c ]) = J(£,,°", [E c ])/J— x(°~~ 
£; [E*]), and J(£; [E c ]) is replaced with J(cr — £; [££]); here, £7* symbolizes the reversed 
profile, E*(£) = E c (a — £). The function Xi(£> [-^c]) is discontinuous at the right contact 
£ = a, and the contact resistance is located there. Since the quasi- Fermi level Ep(£) and the 
thermo-ballistic current «/(£) depend on the chosen location of the contact resistance, they 
are to be considered as auxiliary functions that have no immediate physical meaning. They 
are, nevertheless, of great help in the interpretation of the results, and will be considered 
in this light in the discussion of the examples of Sec. V. On the other hand, the reduced 
resistance x, which is a physical quantity entering the current-voltage characteristic, must 
be unique. This requirement can be met in a natural way by symmetrization: since we have 
Xi(0; [E c ]) = x(°"i 5 we take x m Eq. (|55j) in the symmetrized form 

X = I Me; [E c ]) + X i(cr; [E c ])} = \ { X (a; [E c ]) + X (a; [E*])} , (58) 

which replaces Eq. ()34j) . For constant profile, in which case E* = E c , Eq. (|55j) reduces to 
Eq. In analogy to Eq. (|4"7|) . it is convenient to write 

X = 1 + ko- , (59) 

where 

k=±{k(*;[E c })+k(ct;[E* c })}. (60) 

Since a = S/l and the sample length S is included in the definition of the profile E c , the 
quantity x (and thus also k) is solely a function of I and a (symmetrized) functional of E c . 
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D. The degenerate case for zero bias 



In the degenerate simple treatment is possible only in the limit of zero bias. In 

Eq. (J2HJ), we must replace the current </(£',£") of Eq. with 

J d (e', C) = v e N c (3 ln(l + C -W(C^)-^(PM) [E F (0 - £ F (£")] , (61) 

where Ef(£') and E-p(£") deviate infinitesimally from £p(0). Introducing, in analogy to 
Eq. (J22J, the resistance function 

X d (0 = jJ d (0,0, (62) 

we write </ d (£',£") in the form 

</ d (r, n = Jh d (e, ewe) - i\e , o/coi , (63) 

where the profile function 7 d (£', £") * s defined as 

7 d (£' f") = 1+6 I (64) 

We then find for x d (£) an integral equation which is identical to that for x(0 [ c ^- Eqs. (|38|h 
(JMJ), and (g21)], except that the function y(£ , ,Z") is to be replaced with 7 d (<f,£"). The 
zero-bias conductance per unit area, G = \j/V\v->o, is obtained as 

G = e 2 v e N c (3 ln(l + e"^) ^ , (65) 

where 

X d = \{x d (cr;[E c })+x d ^[E:])} (66) 

by analogy with Eq. (|55Jl . 



V. ILLUSTRATIVE EXAMPLES 

In order to illustrate the formalism developed in the previous sections, we now consider 
a number of simple examples. We present results of calculations for the central quantity 
appearing in the current-voltage characteristic, the reduced resistance, as well as for the 
quasi- Fermi level and the thermo-ballistic current. In particular, we study the dependence 
of these quantities on the ratio of mean free path to sample length and on the shape of 
the band edge profile. To assess the effect of the dimensionality of the electron motion, 
we compare results for one-dimensional and three-dimensional transport; considering the 
properties of the probability functions p n (Q [cf. Sec. II], we may expect the results for 
two-dimensional transport to fall in between those for one- and three dimensional motion 
throughout. 

We first examine the case of constant band edge profile. As an example for non-constant 
profiles, we then consider profiles containing a single maximum; this case includes profiles 
that rise or fall monotonically over the length of the sample. The profile function 7(£ / ,£") ) 
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which is given [cf. Eq. (|30j)] in terms of the maximum £") of the profile in the interval 

[£',£"], can be reduced in all these cases to an explicit function of the profile E c (£) itself. 
This in turn allows the functions /(£) and /C(£, £') to be written in a more explicit, simplified 
form. For profiles more complicated than the ones considered here, the profile function in 
general must be determined numerically. 



A. Constant band edge profile 

In the case of constant conduction band edge profile ("flat band edge"), we have 
7(£'>£") = 1- For the coefficient function / and the kernel /C, we then find from Eqs. 
and flUJ), using Eqs. © and ©, 



/(O=Pd(0)-po(O (67) 

and 

fc(e,o=pi(0-pi(e-o. (es) 

respectively. As the functions / and K. are universal functions of the scaled variables £ and 
£', the solution x(0 °f Eq. (f4~2"j) will also be solely a function of £. 

For one- dimensional motion, we find from Eqs. ([67|) and (|68p. using Eq. (JTUJ), 

/(£) = 1 - (69) 

and 

/C(£,0=^'- e -^'), (70) 
respectively. The solution of Eq. (J42|) is then seen to be 

X (£) = l + £/2, (71) 

which is valid in the range < £ < a [cf. Eq. (}4*3*j) ]. and the function «(£) introduced in Eq. 
f!47|) assumes the constant value 1/2. 

Specializing now Eqs. ([49| ) -(foT ]) to the one-dimensional case, we obtain for the current- 
voltage characteristic 

j = -ev c N c e -fiM0)-E*m _J__ (i _ e -^V) ( 72 ) 

1 + 0-/2 



(with ev e N c = 2e//3h), for the quasi- Fermi level 

1 -l-t lo . _. 1 

(73) 



£ F (£) = £ F (0) + - In 



1 _ i±iZi (i _ e -^) 



and for the thermo-ballistic current 

m 



l + Ue^-e-^) . (74) 



J 

For three-dimensional motion, we have from Eqs. (|57)l and dSEJ), using Eqs. (|T7jl and (fTH|). 

/(£) = | - 2£ 4 (£) (75) 
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and 

ic^o = nm')-m-o}, (76) 

respectively. Then, solving Eq. numerically for the resistance function x(0-> we obtain 
the reduced resistance % = x( a )i which determines the current- voltage characteristic via 
Eq. (}4T?j) . The quasi- Fermi level and the thermo-ballistic current can be calculated from 
Eqs. (JHIIJ) and (J5T]) . respectively. 

In order to characterize the behavior of the reduced resistance % for flat band edge profile, 
we show in Fig. 0] the quantity H [cf. Eq. (J60j) ]. plotted against l/S. For three-dimensional 
transport, the curve for H varies smoothly between the limits 1/5 C 1 (diffusive limit), 
where k — > 3/4, and l/S 1 (ballistic limit), where k — > 1; accordingly, there is a smooth 
transition of the current-voltage characteristic from diffusive to ballistic behavior. For one- 
dimensional transport, we simply have ~k = 1/2 [cf. Eq. 1)71))]. The reduced resistance 
is higher for three-dimensional than for one- dimensional transport, owing to the smaller 
effective mean free path in the former case. In fact, in the diffusive limit, we have x ~^ kS/1, 
and the ratio of the reduced resistances for three-dimensional and one- dimensional transport 
is equal to 3/2, which is the inverse of the ratio of the corresponding effective mean free 
paths (cf. Sec. II). This result holds for any profile [cf. Eq. (jA.17j) ]. On the other hand, in 
the ballistic limit the ratio of the quantities k for three-dimensional and one-dimensional 
transport is equal to 2, not 3/2. However, this has only little effect in the reduced resistance 
because the contribution of k in Eq. ()59|) is suppressed as the ballistic limit is approached. 

Also shown in Fig. 0] is a curve for k that has been derived from the results for flat 
profile and three-dimensional transport obtained in Ref.— and that, while coinciding with 
our results in the diffusive and ballistic limits, exhibits small deviations from our curve in 
the range of intermediate l/S. We return to a discussion of this feature in Sec. VI. 

To demonstrate the behavior of the quasi- Fermi level Ep(C,), we show in Fig.|S]the normal- 
ized internal voltage difference V(£)/V [cf. Eq. (|53|)]. plotted versus the reduced coordinate 
x/S = l£,/S. Calculated with (3eV = eV/k B T = 10, the results of Fig. are typical of 
the strong-bias regime. The curve corresponding to the diffusive limit l/S <C 1 reflects the 
expression 

E F (0 = E F (0) + - In 

that follows by combining Eq. (jA.17|) [with S = S] with the integrated form of Eq. (jA.18|) 
and that is independent of the dimensionality of the electron motion. For l/S > 0, the 
curves have a discontinuity at x = [recall that E F (0 + ) ^ Ep(0), and hence V^(0 + ) ^ V(0)], 
which increases roughly proportional to \n(l/S) up to l/S « 1000. This voltage gap is 
associated with the contact resistance at £ = 0. The latter amounts to the total resistance 
in the ballistic limit, but its contribution to the voltage difference V(£) becomes smaller and 
smaller as the diffusive limit is approached. The results for three-dimensional transport are 
below those for one- dimensional transport, for the same reasons as above. For values of (3eV 
smaller than the chosen value of 10, all curves in Fig. |S] would be shifted upwards. 

The thermo-ballistic current J(£) for flat band edge profile, normalized to the physical 
current J, is shown in Fig. @as a function of the reduced coordinate x/S = l£/S. The 
qualitative behavior of the curves for small x/S can be understood by reverting to Eq. (|52j): 
with x(0 + ) = 1 [cf. Eq. (|43))] and x'{0 > (as can be inferred from the universal curve for 
Co = in Fig. 0J), we have </(£)/ ^ > 1 in the vicinity of x = 0. For x in the vicinity of S, we 
cannot argue in this way; however, the drop of J(£)/J below unity in this range is clearly 
in agreement with Eq. (|2~7I) . 



i - 1 (i - e-n 



(77) 
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For one-dimensional transport [cf. Eq. (J7l*|) ]. the function [J(£)/J] — 1 is strictly an- 
tisymmetric with respect to an inversion of the reduced coordinate at x/S = 1/2. This 
symmetry is (weakly) violated in the three-dimensional case. In conformance with the dis- 
cussion of Sec. IV. C, the thermo-ballistic current in the ballistic limit, l/S ^> 1, tends to 
become equal to the physical current over the whole sample, J(£)/J — > 1, independent of 
the dimensionality of the electron motion. The same holds in the diffusive limit l/S <C 1, 
except for the vicinity of the left and right ends of the sample, where J(x — ► 0)/J = 3/2 and 
J {x/S — > 1)/J = 1/2, respectively, for one-dimensional transport, and J(x — > 0)/J = 1.5485 
and J {x/S — > 1)/J= 1/2 for three-dimensional transport. 

This behavior of the thermo-ballistic current is in accordance with the comments made 
in connection with the internal voltage difference. In the ballistic limit, the total resistance 
is equal to the contact resistance, and thermo-ballistic, ballistic, and physical current all 
coincide. In the diffusive limit, the effect of the contact resistance is suppressed, and the 
thermo-ballistic current becomes entirely diffusive and equal to the physical current. How- 
ever, the diffusive component is defined only inside the sample, and the discontinuities of 
</(£) at the contacts for l/S — > reflect the residual influence of the contact resistance. 

B. Non-constant band edge profile 

For band edge profiles E c (£) that are not constant over the length of the sample, the 
reduced resistance x(7, [E c ]) is to be calculated from the general form of Eq. ()42|1. with the 
coefficient function /(£) = /(£; [E c ]) and the kernel /C (£,£') = /C (£,£'; [E c ]) depending on 
E c (£) via the profile function 7(£', £")■ 

Examples with non-constant profiles allow us to study the interplay of the effects of the 
profile structure and of the magnitude of the mean free path. We investigate the behavior of 
the internal voltage difference, of the thermo-ballistic current, and of the reduced resistance. 
In the cur rent- volt age characteristic (J35|) . of course, the main effect of the profile generally 
is the barrier effect represented by the factor e~@ Ep . 

Here, we consider the case of profiles exhibiting a single maximum at Si in the interval 
[0, S], i.e., at £ = a\ — Si/l. From Eq. (J3DJ), we find for the corresponding profile function 



7 « ( ^ n = py( max (£' ; f «))e (ori _ + 7(^)6^' - a^/m")}®^ - n 

+ [7(ai)e( ffl -O + 7We',O)e(e'-^i)]/[W0(e''-^, (78) 



where we have defined 




where 
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For the kernel /C* 1 ) (£,£') corresponding to 7^(00); we get from Eq. 

/C«(^,0 = tfo&oefa - + £ (>) (£,Oe(£ - *i) , (83) 

with 

/c (<) (£,0 = pi(0-pi(o) 



+ [1/7(01 



4) 



(e - Oj*(* - 07(0 + / w - Oft»(^ - 07(0 

J£' 

and 

/c (>) (e,0 = {[(n -Oftfa -O +Pi(n-0 -Pi(£-0][7(*i)/7(0] 

+ pi(0 - pi(o) + £ #"(£" - Ops(0 - 0[7(0)/7(0l} ©(^i - O 
+ |bi(o) -pi(C - 0] [7(0/7(^1)] - (0 - ^(0 - 

+ pi(o - mo - *o - £ <*o(0 - oxo - om")/^)]} ©(o - ^) . 

(85) 



In the limits o\ — > cr and <7i — > 0, we obtain from Eqs. (|82|) and (jHSj) the coefficient function 
and kernel for monotonically rising profile and falling profile, respectively. 

Considering specifically a piecewise linear profile with maximum E® at Si, we write 



where 



so that 



E c (£) = E° c h(lO = E° c h(x) , (86) 



h(x) = i (s _ _ Si) ] x ~ ^ ' (87) 



7(0 = e- com) , (88) 

with c = /3£ c °. 

In the limit o\ — > cr, i.e., for a profile rising monotonically over the length of the sample, 
the coefficient function and kernel entering Eq. (|42|) are given by the function /«) of Eq. (jUJ) 
and by the kernel K^' of Eq. ()84|). respectively. For this case, by numerically solving 
Eq. (|42jl for three-dimensional transport, we have obtained the curves for the normalized 
voltage difference V(£)/V [i.e., for the quasi-Fermi level £"f(0] shown in Fig. [7J and for the 
normalized thermo-ballistic current J(£)/J shown in Fig. [HJ The flat-profile results (cq = 0) 
for three-dimensional transport (cf. Figs. and EJ) are included for comparison. 

Two features of the curves shown in Fig. are notable, (i) With increasing value of the 
parameter Cq, the quasi-Fermi level -Ef(0 becomes less sensitive to a change in the ratio l/S; 
for Co = 10, the curves for the different values of l/S virtually coincide, (ii) With increasing 
value of Co, the curves tend to stay closer to the abscissa, and the sharp rise in the curves is 
shifted progressively to the vicinity of x/S — 1. Both features can be explained by looking 
at expression (|46p. in which the exponential term depending on the band edge profile E c (£) 
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takes the form exp(— co[l — x/S]) for the linearly rising profile [cf. Eqs. and (|57j) ]. For 
Co 3> 1, this term is sizeable only in the immediate vicinity of x/S = 1, where x(Q ~ x(5/Z), 
independent of S/l. Thus, it emerges that in the presence of a pronounced maximum in the 
profile, the entire internal voltage change occurs in the vicinity of this maximum. 

The localization of the voltage change near the maximum of the profile appears to imply 
that the thermo-ballistic current J(£) approaches the physical current. This follows from 
the results shown in Fig. |S] where the normalized thermo-ballistic current </(£)/ J tends to 
approach unity when the parameter Co is increased, regardless of the value chosen for l/S. 

The effect of a variation of the position Si of the maximum is illustrated in Fig. H 
From the numerical solution of Eq. (|42|) with the coefficient function /W of Eq. (jSUj) and 
the kernel KF> of Eq. (E3J) for one-dimensional transport, we have obtained the curves for 
the quantity k [cf. Eqs. (joTJj) and (pOJ)] as a function of Si/S, calculated for c = 1 and 
10 and for various values of l/S (as the curves are symmetric about Si/ S = 0.5, only 
the range < Si/S < 0.5 is shown). Generally, the dependence of k on Si and on l/S 
is weak. For l/S < 1 (diffusive limit), we have K -> (1/2) (5/5) = (1/2)(1 - e _co )/cb, 
independent of 5i [cf. Eqs. ()A.15|) and ()A.16|) ]; the latter is a consequence of the fact that 5 
is determined essentially by the area below the profile, which is independent of the position of 
the maximum. For l/S ^> 1 (ballistic limit), H assumes its largest values when the maximum 
of the profile lies at either end of the sample. Remarkably, k decreases with increasing value 
of the parameter c , i.e., the reduced resistance % ls lowered as the height of the profile 
maximum is increased. However, as mentioned above, in the current-voltage characteristic 
this lowering of the resistance is generally over- compensated by the barrier factor. 

In Fig. El the normalized internal voltage difference V(£)/V is shown as a function of 
the reduced coordinate x/S = l£/S for l/S = 1, Co = 10, and various values of S\/S. The 
value chosen for Co corresponds to a high profile maximum. One observes an abrupt change 
of the voltage difference in the vicinity of the maximum (cf. Fig. 0). The contact resistance 
seems to have been shifted into the sample unto whatever the position of the maximum is. 

VI. COMPARATIVE DISCUSSION 

In order to put the present description of carrier transport into proper perspective, we 
compare it to previous studies attempting to unify ballistic and diffusive transport. 

An early synthesis of the ballistic and diffusive mechanisms for transport across a single 
potential barrier has been put forward in Ref.— (cf. also Ref.—). The barrier maximum is 
enclosed in a narrow ballistic interval, and the transport in the remainder of the sample 
is regarded as diffusive. For this special shape of profile, this ad hoc model gives results 
qualitatively equal to ours. It was first applied to transport across metal-semiconductor 
contacts, and later (cf., e.g., Ref.—) to transport across grain boundaries in polycrystalline 
semiconducting materials. 

In Ref.—, the transition from the ballistic to the diffusive transport regime in high- 
mobility wires has been studied through a transmission approach based on the semiclassical 
limit of the Landauer transport formula- The zero-bias conductance has been obtained for 
homogeneous samples (i.e., for flat potential profile) at zero temperature for two-dimensional 
and three-dimensional transport. The transmission problem has been solved in analogy to 
the treatment of light scattering in a diffuse medium^. Isotropic impurity scattering of the 
electrons has been assumed, and the average transmission probability has been calculated 
within a two-step formalism involving the solution of a Fredholm-type integral equation 
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(Milne's equation) and of a differential equation for the relevant transmission probabilities. 
The approach of Ref.— thus differs from the present one in that it is restricted to the 
case of flat profile at zero temperature, and is based on a different physical picture. As 
a consequence, the results of Ref.— generally deviate from those in the present work; the 
conductance turns out to be the same for one-dimensional transport, but is different for 
three-dimensional transport (cf. Fig. 0J). An interpolation formula proposed in Ref.— has 
been used in Ref.— to calculate the conductance across a single maximum in the band edge 
profile (arising from a grain boundary in a polycrystalline semiconductor). 

In a previous paper by the present authors on the unification of (one-dimensional) ballis- 
tic and diffusive carrier transport, a closed-form expression for the reduced resistance was 
obtained, using ideas similar to those which have led to the concept of the thermo-ballistic 
current in the present paper. In order to achieve this, it was assumed that the current anal- 
ogous to the thermo-ballistic current is conserved at the equilibration points. For arbitrary 
band edge profile, the reduced resistance is then obtained in terms of explicit integrals over 
the profile, so that the interplay of the effects of the profile and the mean free path can be 
studied in a transparent manner. The method of Ref.— has proven useful in the analysis of 
Hall mobility data for polycrystalline silico n 24 i 25 ; for a fully quantitative analysis, however, 
it should be replaced with the present method. 

The main feature of the present description is the partitioning of the physical current 
into the thermo-ballistic current and a background current. This allows us to take into 
account the coupling to the thermal reservoirs, invoking only global rather than local current 
conservation (the use of local conservation of the thermo-ballistic current would actually lead 
to inconsistent results). In Refe i i 13 ' ' , local current conservation was explicitly assumed, 
but it also underlies implicitly the approach of Ref.—. 

In the calculation of the thermo-ballistic current via the solution of the integral equation 
for the resistance function, as it is done in the present work, the contact resistance appears 
in a natural way as an integral part of the total resistance. The resistance function, as 
well as the quasi- Fermi level and the thermo-ballistic current determined by it, are auxiliary 
functions which depend on the choice of reference contact. This is a well-known feature of the 
transmission approach to the conductance problem- In the present approach, the physical 
results are independent of the above-mentioned choice as a consequence of an appropriate 
symmetrization. 

VII. SUMMARY AND OUTLOOK 

We have developed a unified semiclassical description of ballistic and diffusive carrier 
transport in parallel-plane semiconductor structures. A novel theoretical formulation has 
been achieved by introducing the concept of a thermo-ballistic current, which combines in 
a consistent manner the ballistic motion of carriers in the band edge potential and their 
thermalization at randomly distributed equilibration points. Global conservation of the 
thermo-ballistic current has been established by identifying its average over the sample 
length with the physical current, for which the current- volt age characteristic has been derived 
in terms of a reduced resistance. The formalism is applicable to the treatment of transport 
across arbitrarily shaped band edge profiles. 

Calculations of the reduced resistance have been carried out for simple band edge profiles. 
It has been found that generally the transport mechanism changes smoothly from ballistic to 
diffusive as the value of the ratio of mean free path and sample length changes from large to 
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small. The introduction of a pronounced profile structure appears to weaken the dependence 
on this ratio. The dimensionality of the electron motion turns out to have fairly small effect, 
giving rise to a value of no more than 3/2 for the ratio of the reduced resistances for three- 
dimensional and one- dimensional transport. By using, in one-dimensional calculations, the 
effective mean free path of the three-dimensional case, one may closely approximate the 
results of three-dimensional calculations. Analyses of the local quasi-Fermi level and of the 
thermo-ballistic current have provided detailed insight into the interplay of the ballistic and 
diffusive transport mechanisms. 

The description may be extended in several ways. So far, we have regarded the band 
edge profile as given, so that the integral equation for the resistance function is linear; this 
is justified for weak bias. In the general case, the equation is nonlinear and must be solved 
self-consistently together with the Poisson equation for the band edge profile. Tunneling 
effects can be incorporated in the formulation by changing the classical expression for the 
transmission probabilities to the corresponding WKB expressions. Inhomogeneously doped 
samples like, e.g., p-n junctions, can be treated by introducing a position-dependent mean 
free path. This involves a simple generalization of the expressions for the probabilities of 
occurrence of the ballistic intervals. Another possible extension is to include the effect of 
illumination. To this end, the formalism is immediately adapted to the transport of holes. A 
mechanism for the absorption of light must be introduced, which would eventually lead to the 
determination of the photoconductance. Implementing extensions of this kind in the unified 
description of carrier transport should enable one to apply it in the analysis of experiment 
data, e.g., for polycrystalline semiconducting materials where carrier transport is governed 
by sequences of potential barriers arising from carrier trapping at grain boundaries. 

An extension of the present approach of great topical interest would be to include the 
electron spin as an additional degree of freedom. This would allow us to treat spin-polarized 
transport and, thereby, to make contact with the newly developing field of semiconductor 
spintronic n 26 ! 27 ! 28 . There, the drift-diffusion approach has commonly been use d 29 ' 30 to de- 
scribe the transport of spin-polarized electrons. Only very recently, a synthesis of diffusive 
and ballistic spin transport has been attempted 3 ^. 

In conclusion, we point out that the applicability of our approach is limited mainly by 
two features. First, the assumption of a one-dimensional geometry leaves as principal fields 
of application the transport in parallel-plane heterostructures (including two-dimensional 
wires) as well as in polycrystalline semiconducting materials with a virtually one-dimensional 
grain-boundary structure. Second, quantum coherence effects are beyond the scope of the 
present semiclassical approach. Therefore, limits are set to the treatment of transport in 
nanostructures at low temperatures, where the phase coherence lengths play a dominant 
role. 
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APPENDIX: BALLISTIC AND DIFFUSIVE LIMITS 



For given band edge profile, we consider two limiting situations for the mean free path I, 
the ballistic limit and the diffusive limit. To this end, we rewrite Eqs. flHHj), (JHHJ), and (jj2j) 
in a form exhibiting I explicitly: 

rp PX A ' rp? rp rp? 

f( X ;l) = ± Pl ( X /l)+ ^ ^^_^ p2 (( x _ x ')//) 7 ( x ' )X)) (A.l) 

}C(x,x';l) = — -p 2 (^70 j — P2((a; - x')/l)j(x,x') 

+ I ^-^-^Psdx'-x'yi)^"^') , (A.2) 

j - /(*; X (x; + jf — /C(s, x'; Q /) = . (A.3) 
(i) In the ballistic limit l/S 3> 1, we write 

5 

~ a(x) + 6(x)y , (A.4) 
and find from Eqs. ()A.l|) - (jA.3j) . to second order in S/l, 



s js 
III 



+ p 2 (o) 



5 



-1 + 



s dx 

o 



(l-x/5) 7 (x,5) 



a(5) + 6(5) 



/5\ 2 dx 

+ P 2 (0)( T J ^ -[-a;/5 + (l-a;/5)7(5,x)]a(x) = 



Collecting terms of order S/7, we obtain 



a(S) = 1 



while the terms of order {S/l) 2 yield 



W = p 2 (o) + / y a - [7(5, *) - 7 (x, S)] 



(A.5) 



(A.6) 



(A.7) 



As a result, we arrive at 



x(U^c])«l+P2(0) y |- + jf |. (1 - s/S) e" 



■PE?{x,S) r e /3E c m (0,x) _ e /3£ c m (0,S)l 



(A.8) 



For Z/S 1 — > oo, we have x(Z, [i? c ]) = 1, and Eq. (|33j) goes over into the current- voltage 
characteristic in the purely ballistic regime, 



(A.9) 



21 



in this regime, the current (or the conductance) is independent of the band edge profile E c 
(except for the height of its absolute maximum, -E^O, S), in the interval [0, S]), and of the 
dimensionality of the electron motion. 

(ii) In the diffusive limit l/S 1, we write 



B(x) 



(A.10) 



where B(x) is independent of I to leading order. We then find from Eqs. (|A.lj) - (|A.3)l 

rx fix' r 2 

x-2 Po (0)-f(x,x)B(x) + — 

Jo I Jo 



dx' r r' dx" (x" - x'f , d 7 (x",x') 

-p 3 ((x -x )//)) 



+ 



I I 

* dx" (x" - X'^ 2 



I 



I 



p 3 ((x"-x')/l)) 



d^(x" , x') 
dx 1 



dx' 

B(x ) = . 

(A.ll) 



Closer inspection shows that for given x', depending on whether E c (x') > E c (x") or E c (x') < 
E c (x") when x' ~ x", either only the first or only the second of the two integrals over x" 
contributes, either one yielding 



2p (0) f ic'^M^e 
jo dx' 



We then have 
x - 2p (0) 



d 



e PE ? (0,x) e -pE c (x) B ^_ f fa, e PE?{0,x>) r e -^cM| B( , 

Jo V dx' 1 



(A.12) 



. (A.13) 



Taking the derivative with respect to x, we find 



■y\ d 



1 - 2p (0)e-^W — \e^°^B(x)} = 
dx 



so that 



x(S;t)=x(l,[E$ 



B(S) S([E C }) 



I 2p (0)Z ' 
where we have introduced the effective sample length 

S([E C ])= f S dxe -0[E?(O,S)-E c (x)] 

Jo 

Equation ()A.15|) implies for Eq. (|3*3|) 



(A. 14) 



(A.15) 



(A.16) 



2(0 
S([E C ]) 



1-e 



(A.17) 



Thus, the current- voltage characteristic is seen to depend, in an integral way via S([E C ]), 
on the band edge profile E c (x), and on the dimensionality of the electron motion via the 
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effective mean free path (I) = po(0)l [cf. Eq. ([§}]. Expression (jA.17|) is equal to the result 
obtained for the stationary electron current by integrating the drift-diffusion equation 10 

j — fj,n(x) —Ep(x) , (A.18) 
ax 

using Eq. (j2J); here, /i = 2ev e (3(l) is the electron mobility. 

It appears from Eq. (jA.15|) that the criterion for the diffusive limit, l/S <C 1, should be 
more rigorously taken as l/S *C 1. 
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FIG. 1: Schematic diagram showing a sample of length S enclosed between two plane-parallel 
contacts with associated Fermi levels Ep(0) and Ep(S). 
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FIG. 2: Diagram illustrating the four types of ballistic current entering into the right-hand side of 
Eq. (H|). 
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FIG. 4: The quantity k characterizing the reduced resistance x f° r flat band edge profile, plotted 
as a function of l/S. Solid curve: three-dimensional transport; long-dashed line: one-dimensional 
transport; dotted curve: result derived from the results of Ref^. 
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FIG. 5: The normalized internal voltage difference V(£)/V for one- and three-dimensional transport 
across a flat band edge profile, plotted as a function of the reduced coordinate x/S = l£/S for 
various values oil/S (fleV = eV /k&T = 10). The dotted curve labelled "0" corresponds to l/S <C 1 
and has been calculated from Eq. (|77|l. 
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FIG. 6: The normalized thermo-ballistic current J(£)/J for one- and three-dimensional transport 
across a flat band edge profile, plotted as a function of the reduced coordinate x/S = l£/S for 
various values of l/S. 
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FIG. 7: The normalized internal voltage difference V(£)/V for three-dimensional transport across 
a linearly rising band edge profile, plotted as a function of the reduced coordinate x/S = l£/S for 
various values ofl/S and of the parameter co (cq = 0: flat profile [cf. Fig. E] ; fieV = eV/k-QT = 10). 
Note that for co = 10, the three curves corresponding to l/S = 0.1, 1, and 10 coincide within the 
accuracy of the drawing. 
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FIG. 8: The normalized thermo-ballistic current J(£)/J for three-dimensional transport across 
a linearly rising band edge profile, plotted as a function of the reduced coordinate x/S = l£/S 
for various values of l/S and of the parameter co (co = 0: flat profile; cf. Fig. EJ). For better 
distinguishability, the curves corresponding to l/S = 1 have been shifted downwards by 0.5 units. 
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FIG. 9: The quantity k characterizing the reduced resistance x, calculated for one-dimensional 
transport across a band edge profile with a single maximum at £ = o\ = Si/ 1, and plotted as a 
function of Si/S for cq = 1 and 10 and for various values of l/S. 
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FIG. 10: The normalized internal voltage difference V(£)/V, calculated for one-dimensional trans- 
port across a band edge profile with a single maximum at £ = o\ = Si/ 1, and plotted as a function 
of the reduced coordinate x/S = l£/S for l/S = 1, cq = 10, and various values of Si/S. 
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